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Abstract: We study numerically the Schwinger-Dyson equations for the coupled system 
of gluon and ghost propagators in the Landau gauge and in the case of pure gauge QCD. 
We show that a dynamical mass for the gluon propagator arises as a solution while the 
ghost propagator develops an enhanced behavior in the infrared regime of QCD. Simple 
analytical expressions are proposed for the propagators, and the mass dependency on the 
A.QCD scale and its perturbative scaling are studied. We discuss the implications of our 
results for the infrared behavior of the coupling constant, which, according to fits for the 
propagators infrared behavior, seems to indicate that as{q'^) a.s q"^ 0. 
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1. Introduction 

Due to the property of asymptotic freedom Quantum Chromodynamics (QCD) has been 
extensively tested in the regime where high energies are transferred between quarks and 
gluons. On the other hand we could say that we have only a qualitative understanding of 
its infrared (IR) properties. Several nonperturbative techniques have been used to study 
the infrared region, and among these are the QCD simulations on the lattice, which are 
showing strong evidences that the gluon propagator is infrared finite and that the 
coupling constant may freeze at low momenta |^. 

A possible infrared finite behavior for the gluon propagator and the running coupling 
constant has also been determined as a solution of the coupled system of Schwinger-Dyson 
equations (SDE) for the gluon and ghost propagators ^, ^. Other signals for the 
freezing of the coupling constant can also be found in the literature 0, and there is 
an accumulation of phenomenological evidences and possible tests for this soft infrared 
behavior of the gluon propagator as well as for the coupling constant p. 

In principle the Schwinger-Dyson equations provide a powerful tool to study the QCD 
infrared behavior, because they comprehend an infinite tower of coupled integral equations 
that contain all the information about the theory. In practice its intricate structure only 
become tractable when we make some approximations and truncations. To illustrate the 
difficulty of this method we notice that only in the nineties, Curtis and Pennington pointed 
out a truncation scheme, for QED, which is gauge independent and also respect the multi- 
plicative renormalizability [p!o| . At the moment there are only indications that it is possible 
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to do the same in QCD |11|, and, unfortunately, we have to go step by step changing and 
improving the approximations in order to unravel the actual behavior, or at least to obtain 
rough solutions that could be compared to other methods like QCD simulations on the 
lattice. 

Recently we solved the SDE for the gluon propagator in the Landau gauge within the 
so called Mandelstam's approximation |12]. This is an interesting example of how delicate 
are these equations. 

The gluon propagator within this approximation was first shown to behave as 
in the infrared |13|, which was appraised as a clear signal of confinement. However this 
result was discarded by simulations on the lattice ||^]. In our previous work jl^ assuming 
a different trilinear gluon vertex and renormalization of the final equation, as suggested 
by Cornwall Q, and still in agreement with the Mandelstam's approach, we obtained an 
infrared finite gluon propagator. Our result showed how the approximations may change 
the solutions of the SDE. Of course, although the result is very instructive it is not complete 
because in this approximation the ghosts are neglected. In this work we will improve this 
approximation with the inclusion of the ghosts fields. 

We will solve the coupled SDE for the gluon and ghost propagators in the Landau 
gauge. Fermions will not be included at this level. As in ref. we follow Cornwall's 
prescription to deal with the equation for the gluon propagator. We obtain a numerical 
solution indicating the generation of a dynamical gluon mass without the introduction of 
any ansatz for the solutions. The integral equation for the gluon propagator is clearly 
consistent with a massive gluon polarization operator and its ultraviolet behavior is also 
consistent with perturbative QCD. The distribution of our paper is the following: In section 
II we build up the system of coupled equations. In section III we discuss the angular 
approximation that we perform in the integral equations, which simplifies considerably 
the amount of numerical work to solve the equations. Section IV contains a discussion of 
the renormalization procedure and in section V we check the ultraviolet behavior of the 
equations in order to compare it with the predictions of perturbative QCD. In section VI 
we present our numerical results and discuss its implications for the infrared behavior of 
the coupling constant. Our last section contains the conclusion. 



2. The coupled SDE for the gluon and ghost propagators 

The SDE for the gluon propagator in pure gauge QCD is shown diagrammatically in figure 
(|l|). The gluon propagator is written in terms of itself, the full 3 and 4-point gluon vertex 
r^jyp and r^jypo-; and also the full ghost propagator and the gluon-ghost coupling. The ghost 
propagator depends on itself, on the gluon propagator and also involves the gluon-ghost 
coupling. 

The truncated renormalized SDE for the gluon propagator, in the Landau gauge and 
at one-loop level, are shown in figure @i in the same approximation of ref. and it can 
be written in the Euclidean space as 
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Figure 1: The complete Schwinger-Dyson equation for the gluon propagator without quarks. The 
spiral lines represent the gluon fields and the dashed lines the ghost fields. The black blobs indicate 
the full gluon propagators while the white ones represent the full vertices. 



(k) 



1 



(27r) 



-ro (k 

4 /j.pay'^^ 



-P, q)Da(s{q)DpaTi3„i,{-q,p, -k), (2.1) 



where p = k+q, D^^^{k) and rj^^p(p, q, k) are the gluon propagator and three gluon vertex at 
tree level, V ^yp{p, q, k) and Gy are the full three gluon and gluon-ghost vertices respectively, 
and we use the color factor C2 = 3. Note that all terms with four-gluon vertices were 
neglected. One of these contributions is the momentum independent tadpole which can be 
eliminated by one appropriate choice of the momentum projector as discussed in second 
paper of ref. [13|. 



The full gluon propagator that enters into eq.(2.1) is expressed by 



D'"'{k) 



fc2 



A;2 
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where 2{k'^) is the dressing of the gluon propagator. When Z{k'^) = 1 we recover the 
perturbative expression of the gluon propagator at tree level. Therefore this function 
measures the transition from the nonperturbative to the perturbative regimes as its value 
changes with the scale. 

It will be useful to introduce here the function D{k'^) that can be expressed in terms 
of the dressing of the gluon propagator, Z{k'^), 
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Figure 2: Diagrams for the coupled gluon-ghost system of Schwinger-Dyson equations. 



Dik') = (2.3) 
The full ghost propagator, Ddk), can be defined as 

Dcik) = (2.4) 

where T{k'^) is the dressing of the ghost propagator. 

The renormalized SDE for the ghost propagator, shown in figure reads 



D^\k) = -Z^k^ + g'C2Zi j ^^ik^D^Uk - q)G,{k, q)DG{q). (2.5) 

Note that in the above equations, we have already introduced Z3, Z3, Zi and Zi which 
are respectively the renormalization constants for the gluon propagator, the ghost propa- 
gator, the three gluon and the gluon-ghost vertices, which are needed in order to render our 
coupled integral equation system finite. Once the full three gluon and gluon-ghost vertices 
are known, we have a closed coupled system for the gluon and ghost propagators that can 
be solved numerically. 

The construction of these vertices is based on their respective Slavnov- Taylor identities 
and its form was discussed in ref. Q. Apart from a group theoretical factor (g/"*'^) we can 
express the full three-point vertex function in the following way 

^^iup{p, q, k) = -A+(/, q^, k'^)S^^i{p - q)p - A^ip^, q^, k^)6^ui{p + q)p 
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"^ "^ -q2^ - Pv%)i{P - q)p + c.p., (2.6) 

where 

The three gluon vertex at tree level can be recovered if we assume that 2{k'^) = 
!F{k'^) = 1 and in this way we obtain A_{p'^ ,k'^) = and g^, A;^) = 1. 

The full gluon-ghost vertex can be written as [Q] 

GM 1) = 'i^y^ + 'P^^ {y^ - 1) • (2.8) 

Introducing into the eq.(|2.5|) the eqs.( |2.2| ), ( |2.4| ) and ( |2.8D , we obtain 

- Z^-g C2Z1 / — - {kV{p)q) , , ,0 o + ^FTT^T " ^ ' (2-9) 



^(A;2) ^ " V (2^)4 ' k^p^q'^ KHq'^) H^^) 

where p = k — q and 

pM^(p) =5'^'^-^, (2.10) 

is the transversal projector. 

Equations ( |2.lD and (|2.9D are the coupled SDE equations for the gluon and ghost 
propagators in the approximation of ref. 0] that we want to solve, following the same 



steps performed in ref. |12| 



3. The angular approximation 

In order to solve the coupled system of SDE, compound by the eqs.( |2.1| ) and (^^), shown 
in the previous section, the first step is to perform the angular integration of the coupled 



system of integral equations. In our previous work |12|, where we considered only the SDE 
for the gluon propagator and used a specific form for the three gluon vertex, it was simple 
to compute analytically the angular integration without any approximation. However, in 
the present case this is not possible anymore and we perform an angular approximation Q 
as we explain in the following. 

When g2 < p we have that Tip"^) = T{{k - qf) J='{k'^) and Z(p2) ^ Zik"^), that 
preserves the correct limit when ^ . On the other hand when we have > k"^ we 
can set J^{p^) ~ ^{k"^) — ^ ^{o^)- According to this approximation the ghost dressing can 
be written as 



Z3 — ^AZi 



(3,1) 



where A = g'^(^^)/167r^ and the momentum integration in the SDE runs from zero to 
infinity, and a ultraviolet cutoff A is introduced to deal with the upper limit of the integrals. 
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The above equation lead us to 



(3.2) 



while the equation for the gluon propagator now reads 
1 




where we used the projector 

n'^'ik) = S"" - 4—^, (3.4) 



It is important to mention that to obtain eq.(3.3), it was neglected one term from the 
three-gluon loop when < fc^, as discussed at length in the second and third papers of 
ref.i. 

The above equations are identical to the ones obtained by the authors of ref. [Q] in 
their first calculation of the coupled SDE, and they can be viewed as a natural extension 
to the gluon SDE in the Mandelstam approximation, since the only difference between the 
three gluon loop contribution in the previous case and the coupled system is that the gluon 
dressing function Z was replaced by the product ZG. 

Our system of equations will differ from the one obtained in ref. Q after the renor- 
malization procedure. One natural question that can arise is how much our solutions will 
depend on the angular approximation introduced in this section. Obviously, this can only 
be answered after improving this angular approximation, what we hope to check in the 
future, however we can say in advance that it is known that in the case of ref. Q this 
approximation did not introduce a qualitative change of the solution, but it does produce 
a quantitative change, specially in the infrared fixed point value of the coupling constant 

As we will discuss in detail in the section V, this angular approximation was built 
in order to recover all the parameters - the gluon and ghost anomalous dimensions and 
the first order coefficient of the Callan-Symanzik (3{g) function - that describe the known 
perturbative behavior of the coupling constant. However we believe that this approxima- 
tion, despite the fact that it is sucessful in describing the perturbative region [Q] , can bring 
some numerical uncertainty in the infrared behavior, which will be reflected in the ratio 
mo/AQCD, which is a parameter, to be introduced in the section VI, that measures the 
value of the gluon propagator at zero momentum. 

4. Renormalization 

As quickly mentioned in the section III, we can link the unrenormalized propagators and 
vertices with renormalized ones introducing the multiplicative renormalization constants 
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in the following way 



DZ.{q\P^) = Z^{^l\^^)D^,{q\^?) 



-\nr ( „2 



(4.1) 



where the superscript nr denotes the nonrenormalized quantity and // is our renormalization 
scale. The renormalization constants, in the Landau gauge, are connected through the 
following relations 

Zi = ZgZ^"^ ; Z\ = ZgZ\l'^Z'>,. (4.2) 



Furthermore we have one more identity, which is satisfied only in the Landau gauge 

To determine the renormalization constants we will impose the following condition on 

the ghost dressing ^(/i^) = 1, where \i is chosen in the high energy region, i.e. \i » Aqqd, 

where Aqcd is the QCD scale. This procedure is explained in our previous work [12| and 

is usual when dealing with the SDE |Q] . 
According to this we obtain 

H^^) = 1, 
where Aj^{x) is given by 



1 - — ^^(^'), 



(4.3) 



1 r^'' dv 

-Z{x)T{x) + / ^Z{y)T{y) 

^ J X y 



(4.4) 



where we substituted fc^ by x and q^ by y. 

We finally obtain the renormalized expression for the ghost SDE 
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(4.5) 



where J-{x) in eq.([4.5|) is now the renormalized ghost dressing. 

In the case of the gluon propagator we can rewrite the renormalized eq. (|3.3| ) in the 
following compact formula 



Z{x 



,-1 



Z3 + \ZiB{x) + \C{x) 



where 



^ dy ( 7y 
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-7+-jz{y)T{y), (4.7) 



and 



C(x) 



dy ( 3y 



X \2x 



Hy)Hx) + 



2y 3 ^ ' 



(4.8) 
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The Z3 renormalization constant role is to eliminate the divergent terms of the gluon 
SDE (eq.( [4.6[ )). Therefore we can add all the potentially divergent terms imposing that 

1 r^^ X r^^ 

xZ-s + - A / dy-T\y) - 7A / dyxD{y)T{y) =R + x (4.9) 

^ Jx y Jx 

where R can be determined through the gluon propagator renormalization condition 

D-\iJ.^) = fi'^ or Z{iJ.^) = l, (4.10) 

where, again, it is worth remembering that the scale /i^ is chosen in the perturbative region. 
It is important to stress at this point that the subtractive renormalization that we have 
done above is the same that has been prescribed by Cornwall many years ago and 
has been explained in detail in ref. 113]. This renormalization allows a massive solution 



for the gluon propagator. Moreover, as explained in |^] and |12], this approach does not 
break the Slavnov- Taylor identity involving the gluon propagator and the trilinear gluon 
vertex as long as we add to the full triple gluon vertex massless pole terms that have been 
usually neglected in these equations. These terms do not modify the ESD but promote the 
consistency with the Slavnov- Taylor identities. Finally note that we have not discarded 
any term from the unrenormalized gluon SDE, but just absorbed the divergent terms in the 
renormalization constant Z3, and that this constant in eq. ([4.9| ) is proportional do 1 plus a 
function of ji^ and A^, which is compatible with the expected weak coupling expansion for 
this renormalization constant. 

In eq.( f4.6D besides the renormalization constant Z3 we have some terms multiplied by 
the constant Zi that comes from the trilinear vertex. We could eliminate this constant in 
a rigorous way using the relations of eq.(|4.2|) such as 



Z'^ 

Zi = ^. (4.11) 

Z3 

which follows from the identity Zi = 1. However it was shown in ref. |Q] that such 
procedure may destroy the perturbative aspects of the solution within this approximation. 
As mentioned in our previous work, ref. [|l2|, the SDE renormalization procedure is really 
an intricate subject and there is not a recipe to deal with these renormalization constants 
in the nonperturbative region. In such a situation the best that can be done is to go step 
by step and analyse the consequences of each choice, and we decided in this work to choose 
the simplest case discussed in ref. [§] , Zi = 1 . Using this relation to determine the constant 
R defined in eq. (|4.g| ), we obtain the following expression 

R=-^lfdy (I^ -'if.-l) yDiy)Hy) - \x[ dyf,H,^)Hy) 

- A / dyl^yD{y)Hy) + (4.12) 



We finally obtain the renormalized SDE equation for the gluon propagator 
Z)-i(x) = i? + X + A - y I - yD{y)ny) + dy^H^)Hy) 
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+ A / dy~yD{y)F{y) - (4.13) 

Jx 3 

where we recall that A = {p?) / IQ-k"^ , x = k"^ and y = . 

It is this last equation, eq.( 4.13|) , together with the one for the ghost propagator, 



eg. ( [4. 5]) , that will be solved numerically. Note that in eq.( [4.13| ) as x — > the inverse 
propagator goes to a constant value, i.e. it shows the presence of a dynamically generated 
mass. 

5. Ultraviolet Behavior 

Before solving numerically the coupled SDE it is important to show that they reproduce 
the QCD perturbative behavior in the high momentum region, confirming the consistency 
of the renormalization procedure. In order to do so we substitute the perturbative behavior 
of the gluon and ghost propagators in their respective SDE and keep only the ultraviolet 
dominant terms in each equation. The solution that comes out are 



7oa(^ ) 1 ( X 



D-\x)=x[l + ^L^\n\-^]] , (5.1) 



and 



«-W=-(l + ^ln(^)). (5.2) 



where 7^ = 13/2 and 6'q = 9/4. 

As happened in the case of ref. |jl^ (see also the discussion of ref. Q) the momentum 
behavior of the gluon and ghost propagators is the one expected from perturbative but in 
addition we recover, in this approximation, the correct perturbative values for the anoma- 
lous dimensions (70 = 13/2 and 5q = 9/4 ). Furthemore, we can also obtain the correct 
value for the first coefficient of the (5 function (/3q). 

Starting from the eqs. (13) and (U) we can determine the running coupling constant 



and /3q. The following relation is satisfied by the renormalized coupling constant 

where the renormalization constants can be obtained from the gluon and ghost dressing 
functions 

z^'{x,K^) = z^{^Ji\^^)z{x,^?) 

r''-{x,A^) = Z,{f,\A')T{x,f,^), (5.4) 
leading to the following expression for the coupling constant 

a,(x,A2) = as{ti^A^)Z{x,i^^)J^\x,fi^). (5.5) 
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Inverting the above definition of tlie coupling constant, using the perturbative defini- 



tion of the renormalized gluon and ghost functions, which are giving by eqs.(5.1) and (| 
divided by x; and keeping only terms 0(a(^^)), we can obtain the following expression 

1 - ^ .(ii^4^y p.o) 



as{x) cts{^^'^) 47r \A*^/ 

where /3q = Jq + 25q = 11, which is identical to the perturbative value (/3o = 11)- Expressing 



the eq.(5.6) in a more familiar way, lead us to 



An 

as{x) = r-, (5.7) 

where Aqcd is the usual QCD scale. 
6. Numerical Solution 

We can now turn to the numerical calculation to solve the gluon-ghost coupled-system 
expressed by the eqs.( |4.5D and ( [4.13| ). The solution comes out when we apply an extension 



of the same iterative numerical method used in ref. |12|. We start defining a logarithmic 
grid for the x and y variables in order to perform the integration from the deep infrared 
region to the high energy momenta. This grid is split into two regions, the infrared one 
[0,/i^], and the other corresponding to the perturbative region [/i^,A^]. This splitting is 
needed to impose the renormalization conditions on the dressing functions at the scale 

We then provide initial guesses for the functions D{x) and T{x) and use them to 
generate the coefficients of the cubic spline interpolation which will produce the values of 
these functions in terms of the argument y. These ones are used in the right hand side of the 
equation for computing the integral through Adaptive Richardson-Romberg extrapolation. 

The initial guesses are compared with the numerical results that were obtained after 
the integration. The convergence criteria to stop running the numerical code is that the 
difference between input and output functions must be smaller than 10~^, otherwise these 
new numerical results will feed again the right-hand side of the SDE equations and re-start 
all the procedure until the convergence criteria is satisfied. Using this method we verified 
that our results are independent of the starting guesses for D{x) and J^{x). Our input 
data are the renormalization point, /u^, and the value of the coupling constant defined at 
this scale, as{^j?) = (7^(^^)/47r, and this will determine the value of Aqqd through the 
eq.( |5.7D . We vary both parameter fj,'^ and as(/i^) to scan how our solutions depend on the 
renormalization point. 

We compute the gluon propagator D{x) for various sets of //^ and as(/Lt^) within the 
range [10 GeV^, 30 GeV^] and [0.22,0.30] respectively as shown in figure (^). We can check 
through the table(|^) that these values correspond to values of Aqcd parameter from 
321 MeV to 815 MeV. The range of momenta chosen for the scale fj? is consistent with 
a perturbative scale and is within the window of squared momenta 10~^ to 10^ GeV^, 
which is the maximum range that our calculation can cover without loss of precision in the 
infrared region. 
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q'[GeV'] 



Figure 3: Gluon propagator, D{q^), as function of momcntmn for different scales. Tfie line 
+ square curve was obtained when Q;s(10GeV^) — 0.25 which corresponds to Aqcd — 321 MeV, 
while in the line + triangle curve, as(30GeV^) — 0.30, which leads to Aqcd = 815 MeV. The 
shadowed area delimits the curves with Aqcd varying within the range [321 MeV, 815 MeV]. The 
central curve (line + circle) was obtained when we fix the renormalization point, /i^, at the bottom 
quark mass, = (4-5)^ GeV^ with the central value as{mf) = 0.22. 
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Figure 4: Ghost propagator, Da{q^), as function of momentum q^ for the same scales shown 
in the figure (|). The fine + square was obtained when Q;s(10GeV^) = 0.25 which corresponds 
Aqcd — 321 MeV, while that in the line + triangle curve, as(30GeV^) = 0.30, which leads to 
Aqcd = 815 MeV. The line + circle was obtained when we fix the renormalization point, at 
bottom quark mass, ml = (4.5)^ GeV^ with the central value as{ml) = 0.22. Such curves are 
obtained varying the Aqcd parameter within the range [321 MeV, 815 MeV]. 

Independently of what is the Aqcd value we can see that all curves in figure (^) develop 
the same perturbative behavior, on the other hand these curves split in the infrared region 
going to different values in the limit when the momentum q"^ goes to zero. 

For the same set of parameters we have also plotted, in the figure (Q), the behavior of 
the ghost propagator as a function of the momentum q'^. As we can see, for the ghosts, 
the infrared behavior has a slighter dependence on the renormalization point value than 
the gluon propagator. In the overall, all curves develop the same behavior from the deep 
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Figure 5: The behavior of the gluon propagator, D{q^), obtained through the numerical solution 
of the Schwinger-Dyson equations, when as{n^) — 0.22 at /i^ = = (4-5)^ GeV^, together with 



its ultraviolet behavior given by eq.(5.1) where Po — 11. 



infrared to the high energy region. 

In the figures showing the behavior of the gluon and ghost propagators, figures (^) and 
(Q), the most representative curve is the one where the renormalization point, is set at 
the bottom quark mass, = (4.5)^ GeV^ with as{mf) = 0.22. In this case we obtain a 
value for Aqcd = 335 MeV. We are going to concentrate only in this curve to analyse the 
perturbative and also non-perturbative regimes, however it is important to keep in mind 
that it does not really matter what curve we choose to study the high energy regime, since, 
as mentioned before, all of them have the same perturbative behavior. 

For the input discussed above, we plot in figure (|5|) the gluon propagator, D{q^), 
together with its ultraviolet behavior described by eq.(^.l|) while, in figure (|6|), we compare 
the ghost propagator with its ultraviolet behavior expressed by eq.( |5.2| ). It is interesting to 
note that, in this latter case, the behavior of the nonperturbative ghost propagator is not so 
much different from its perturbative behavior, since the major difference starts happening 
only for (f' values less than 2 x 10^^ GeV^. 

Despite the fact that the ghost fields are important to warrant the gluon transversality 
in the perturbative region and also to recover the first order beta function coefficient and 
anomalous dimension of the propagators, the above result suggests that neglecting this 
field, as happens in the Mandelstam approximation, based on an extrapolation from its 
known small contribution in the perturbative region can be a reasonable approximation. 

In order to obtain a value for the gluon propagator near (7^ = 0, we propose a very 
simple fit for our numerical data, which in Euclidean space can be written as 

and where the dynamical gluon mass M^{q^) is described by 

M\e) = (6.2) 
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Figure 6: The behavior of the ghost propagator, DQ{q^), obtained through the numerical solution 



of the Schwinger-Dyson equation, when as(/i^) — 0.22 at ^j? 
its ultraviolet behavior given by eq.(|5.2|) where — H- 
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Figure 7: Numerical solution for the gluon propagator, D{q^), versus momentum for as{m1) 



0.22. We compare this numerical solution with the fit given by eq.(O), where TOq = 0.99 GeV^ 



The reason for this fit was discussed in our previous work |f2|. In figure (Q) the nu- 
merical solution for the gluon propagator, determined with the renormalization point fixed 
at the bottom quark mass, is quite well adjusted with the above fit for mg = 0.99 GeV^. 
In table (2) the values of ttiq, utilized in the eq.( |6.2D , are shown for each value of the renor- 
malization scale and coupling constant. The value of ttt-q itself is not important because 
it is linked to the scale Aqcd through the running of the coupling constant. What really 
matters is the analysis of the ratio ijiq/Aqcd which ranges from 1.72 to 2.97 and in agree- 
ment with previous determinations for this ratio It is clear that different fits will 
give slightly different values for the gluon mass, but other choices, as long as they respect 
the correct asymptotic limits, do not differ appreciably from the values we quoted above. 
Furthermore the ratio is also relatively stable against variations of the renormalization 
point. 
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{(3 = 11) 






0.22 


20.25 GeV^ 


335 MeV 


0.99 GeV^ 


2.97 


u.zo 


1 n 

iU Kjid V 


ozi ivie V 


U.OO yjQ V 


Z.ol 


0.25 


20 GeV^ 


455 MeV 


1.10 GeV^ 


2.31 


0.25 


30 GeV^ 


557 MeV 


1.65 GeV^ 


2.31 


0.30 


10 GeV^ 


471 MeV 


0.65 GeV^ 


1.72 


0.30 


20 GeV2 


666 MeV 


1.31 GeV2 


1.72 


0.30 


30 GeV^ 


815 MeV 


1.96 GeV^ 


1.72 



Table 1: Values of the renormalization point, /i^, and coupling constant, a(/i^), used as input data 
in the eqs.( |4.5| ) and (4.13). These values fix automatically, through eq.(5.7), the Aqcd scale. In the 
third column, we have the values of Aqcd computed with our perturbative value of /3o — H- The 
values for the ratio mo /Aqcd are also shown in the last column where the uiq values are defined 
through eq. (|6.2| ). 



We can also discuss the behavior of the coupling constant in the infrared. In order to do 
so we start remembering that the running coupling constant is given by eq. (|5.5| ), where the 
gluon and ghost dressing functions are utilized as input. Note that if we assume the fit given 
by eq.(6.1) for the gluon propagator and roughly approximate the ghost dressing function 
by 1, because in the full region of momenta its behavior is almost the perturbative one as 
can be seen in figure @, we clearly obtain a vanishing coupling in the deep infrared region. 
Since the major difference between the nonperturbative and the perturbative behavior of 
the ghost propagator starts happening for values less than 2 x 10~^, instead of the 
rough approximation J^{x) ~ 1 in the low momenta region we introduce the following fits: 
J^{x) = (x/mg)'' for the ghost dressing function, and Z{x) = {x/{x + m^))^ for the gluon 
dressing function, in the region delimited by < 0.1 GeV^. These exponents provide 
additional freedom to our fits, and we are able to fit the data with r] = —0.04 and 6 = 0.98. 
These figures confirm the same vanishing behavior of the coupling constant in the deep 
infrared region previously examined, in the context of lattice QCD simulation, by the 
authors of Ref . IT^l . 



7. Conclusions 



We solved the coupled system of Schwinger-Dyson equations for the gluon and ghost prop- 
agators in the Landau gauge in the case of pure gauge QCD. The new ingredient in our 
approach is that we use a renormalization prescription formulated by Cornwall many years 
ago which allows for a dynamically generated mass. As explained in a previous work 
this solution is compatible with the Slavnov- Taylor identities when a new piece containing 
massless poles is added to the triple gluon full vertex. The renormalization constant has 
the form Z3 = 1 + /(/i^, A-^) which is consistent with the weak coupling expansion of this 
constant. The ratio between the dynamical gluon mass and Aqcd is also consistent with 
previous determinations. 
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To obtain the SDE solutions we do not need any ansatz for the asymptotic equation 
and we are able to solve them numerically in a quite large range of momenta. The numerical 
solutions are quite stable and we could, as we did in ref. [^], compute an effective potential 
for composite operators to show that we also obtain a reasonable value for the gluon 
condensate, although this calculation is a little bit redundant because the gluon propagator 
do not differ appreciably from the one that we obtained in the previous work. It seems that 
in the present case the inclusion of ghosts induces only a minor numerical modification of 
the previous result. We present simple fits for the gluon and ghost propagator and discuss 
the infrared behavior of the running coupling constant. These fits allow us to study the 
behavior of the running coupling in the deep infrared region, and indicate that the running 
coupling goes to zero when — > 0. 

There are many points that still must be improved in the present approach. We need 
to compute the equations without the angular approximation. A different approximation, 
other than Zi = 1 must be tested. The inclusion of fermions and the tadpole diagram are 
among the many other modifications that must be considered in the future. 
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